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Abstract. We present the results of a new anal- 
ysis of data taken in 1998-2002 for a search for 
high-energy extraterrestrial neutrinos. The analysis 
is based on a full reconstruction of high-energy 
cascade parameters: vertex coordinates, energy and 
arrival direction. Upper limits on the diffuse fluxes 
of all neutrino flavors, predicted by several models 
of AGN-like neutrino sources are derived. For an 
E -2 behavior of the neutrino spectrum, our limit is 
E 2 F l/ (E) < 2.9 x 10 7 cm 2 s 1 sr 1 GeV over a 
neutrino energy range 2 x 10 4 ^ 2 x 10 7 GeV. This 
limit is by a factor of 2.8 more stringent than a limit 
obtained with a previous analysis. 

Keywords: high-energy neutrinos, neutrino tele- 
scopes, Baikal 

I. Introduction 

High-energy neutrinos are likely produced in violent 
processes in the Universe. Many theoretical models pre- 
dict that neutrinos are generated by hadronic processes 
within high energy astrophysical sources such as active 
galactic nuclei (AGN), supernova remnants or gamma 
ray bursts. Individual sources might be too weak to 
produce an unambiguous directional signal, however the 
total neutrino flux from all sources could produce a 
detectable diffuse neutrino signal. Astrophysical neutri- 
nos generated in top-down models are, by definition, 
of diffuse nature. To date the highest sensitivities to 
diffuse neutrino fluxes in a range 10 TeV -j- 100 PeV 
are achieved with the NT200(Baikal) fl] and AMANDA 
El, neutrino telescopes. 



The Baikal Neutrino Telescope NT200 is operating in 
Lake Baikal at a depth of 1.1 km and is taking data 
since 1998. Since 2005, the upgraded 10-Mton scale 
detector NT200+ is in operation. Detector configuration 
and performance have been described elsewhere 0J, 
0, 0, Q. Due to high water transparency and low 
light scattering, the detection volume of NT200 for high 
energy v e , and v T events significantly exceeds the 
instrumented volume. Our previous analysis |1 1 of 1038 
live-days data, collected in the years 1998-2002 with 
NT200, has allowed to set the limits on diffuse neutrino 
fluxes predicted by several theoretical models. Here we 
discuss a new analysis method which is based on an 
energy and space-angular reconstruction of high-energy 
cascades and present limits on diffuse neutrino fluxes 
which are improved by a factor of about three over the 
previous ones. 

II. The analysis method 

The BAIKAL survey for high energy neutrinos 
searches for bright cascades produced at the neutrino 
interaction vertex in a large volume around and below 
the telescope. The main background source are atmo- 
spheric muons, with a flux 10 6 times higher than that 
of atmospheric neutrinos. We select events with high 
multiplicity of hit channels N^, corresponding to bright 
cascades. To separate high-energy neutrino events from 
background events, a cut to select events with upward 
moving light signals has been developed. We define 
for each event t m [ n = min(^ — tj), where ti, tj are 
the arrival times at channels i,j on each string, and 
the minimum over all strings is calculated. Positive 
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TABLE I 

Expected number of events N model , detection energy range which contains the central 90% of the expected signal 

AE 90 %, MEDIAN ENERGY OF THE EXPECTED SIGNAL E v , AND MODEL REIECTION FACTORS T] = n 90 %/N m odel FOR MODELS OF 

ASTROPHYSICAL NEUTRINO SOURCES. 





BAIKAL 


AMANDA 


Model 


N model{Ve +V^,+V T ) 






n 90%/ N mo del 


n 90%/ Nmodel 


S05 


0.7 


100 TeV -r 30 PeV 


2 PeV 


3.4 


1.6 


Pp7 


4.4 


320 TeV -r 160 PeV 


6 PeV 


0.5 


0.3 


M pp + P7 


1.7 


20 TeV -r 500 PeV 


15 PeV 


1.4 


1.2 


MPR 


1.4 


160 TeV -r 100 PeV 


3 PeV 


1.8 


0.9 


SeSi 


2.4 


1 PeV -f 50 PeV 


10 PeV 


1.0 





and negative values of t m ; n correspond to upward and 
downward propagation of light, respectively. We require 



Vn > -10ns. 



(1) 



This cut accepts only time patterns corresponding to up- 
ward traveling light signals. It rejects most events from 
brem-cascades produced by downward going muons 
since the majority of muons is close to the vertical; they 
would cross the detector or pass nearby and generate 
a downward time pattern. Only few muons with large 
zenith angles may escape this cut and illuminate the 
array by their own Cherenkov radiation or that from 
bright cascades from below. 

The energy spectrum of neutrinos from galactic and 
cosmological sources or from the decay of topologi- 
cal defects is expected to have a significantly flatter 
shape than the spectrum of atmospheric muons and 
neutrinos. This gives rise to different Nhit and cascades 
energy distributions. Results of a search for high-energy 
neutrinos, based on Nhit as a rough indicator of the 
energy deposited in the effective detection volume were 
published in (Tj. Here we present results of an extended 
analysis which is based on a full reconstruction of 
cascades parameters. 

A. Cascade parameters reconstruction 

We use a two step procedure for cascade parameters 
reconstruction which is applied to events with > 5 hit 
channels. At the first step, using the time information of 
hit channels, the coordinates of the cascade vertex are 
reconstructed by minimizing 



x 2 t 



1 



N hi 



t (Ti(x,y,z,t ) - Uf 



(N hit - 4) ^ 



(2) 



Here ti is the time measured by the i-th channel, Tj is 
the expected arrival time of Cherenkov photons induced 
by a cascade with f s h(x,y,z) coordinates, and a ti are 
the timing errors. At the next step, taking into account 
the found coordinates r s h(x,y, z), we reconstruct the 
cascades energy E s h, as well as zenith and azimuth 
angles 8 and tp of cascade axis by maximizing the 
likelihood function 

N hit 

L A = Y[p l (A l ,E s h,n s h{0^))- (3) 

i=l 

Here pi{Ai,E s h,£t s h{Oi<p)) is the probability to detect 
a signal of amplitude Ai by the i-th hit channel. These 



probability functions are calculated with taking into 
account the absorption and scattering of light in water 
and the relative orientations of optical modules and 
cascades. 

B. High-energy neutrino simulation 

The number of expected events during observation 
time T is 

N V = T J dVL j dE sh V eff (ti,E sh )J2 J N APH2 o 
x^$ u (n,E u ,X)dE u , (4) 

X(&) = / Pearth{l)dl, 
JO 

where E„, X) is the flux of high-energy neutrinos 

with energy E v in the vicinity of the detector, ft - 
the neutrino direction, X(U) - the thickness of matter 
encountered by the neutrino on its passage through 
the Earth, E s h - the energy of secondary cascades, 
V e ff(Cl, E s h) - the detection volume. The index Vi 
indicates the neutrino type and k =1,2 corresponds 
to CC- and NC-interactions, respectively. Na is the 
Avogadro number and ph 2 the water density. 

A MC-code is used to solve Eq. with the boundary 
conditions for neutrino fluxes $> Ui (E,0) = A Vi f Vi (E), 
where f Vi (E) is a diffuse AGN-like flux or other pre- 
dicted UHE neutrino fluxes, and a A Vi a normalization 
coefficient. For neutrino interactions we used cross- 
sections from 0. The neutrinos are propagated through 
the Earth assuming the density profile of the Preliminary 
Reference Earth Model [9]. Although a flavor ratio 
k, 1:2:0 is predicted for generic neutrino 



Of Ve.Vp-.l 

fluxes at astrophysical sources, equal fractions of all 
three neutrino flavors are expected at Earth because of 
neutrino oscillations. Throughout this paper we assumed 
a neutrino flavor ratio at Earth of v^.v^.v-y =1:1:1 
and the same shape of energy spectra f v (E) for all 
neutrino flavors, as well as a flux ratio for neutrino and 
antineutrino of v jv = lfl 

The detector response to Cherenkov radiation of high 
energy cascades was simulated taking into account the 
effects of absorption and scattering of light, as well as 

A violation of this assumption (e.g. for neutrino production in jry 
interactions) has a small influence on the result due to the similarity 
of v and v cross-sections in our energy range. 
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Fig. 1. Left: Reconstructed cascade zenith-angule distribution for data (dots) and for MC-generated atmospheric muons (boxes); true MC zenith- 
angule distribution is given as histogram. Right: Reconstructed cascade energy distribution for data (dots) and for MC-generated atmospheric 
muons (boxes); true MC energy distribution is given as histogram. 



light velocity dispersion in water. We also implemented 
the longitudinal development of cascades. For electron 
cascades with E s ^ >2xl0 7 GeV and for hadronic 
cascades with i? sn >10 9 GeV, the increase in cascade 
length due to the LPM effect [ fOj was approximated as 
E 1 / 3 according to ifTTl . 

C. Atmospheric muon simulation 

Downward going atmospheric muons are the most 
important source of background. The simulation chain of 
these muons starts with cosmic ray air shower generation 
using the CORSIKA program 03 with the QGSJET 
lfl3l interaction model and the primary composition and 
spectral slopes for individual elements taken from |[T4ll . 
Atmospheric muons are propagated through the water 
using the MUM program [15|. During passage through 
the detection volume the detector response to Cherenkov 
light from all muon energy loss processes is simulated. 
A total of 1.2 x 10 9 background events, equivalent to 
3671 live days, has been simulated, with standard optical 
parameters of Baikal water. The cascade reconstruction 
procedure was applied to simulated background events 
with Nhit >15 which obey the condition (1). To reject 
events induced by muon bundles the following additional 
cuts were applied: N^u >18, \t <3 and La < 20. 
Zenith-angule and energy distributions of the selected 
events are shown in Fig.l and are discussed in the next 
section. 

III. Results 

Within the 1038 days of the detector live time, 
3.45 x 10 8 events with N^ t > 4 have been recorded. 
For this analysis we used 18384 events with hit channel 
multiplicity N^ t >15 and £ mm >-10 ns. As it was shown 
in 01 the data are consistent with simulated background 
for both t m m and distributions. 



A full cascade reconstruction algorithm (for vertex, 
direction, energy) was applied to the data 11611 . Cuts 
were then placed on this reconstructed cascade energy 
to select neutrino events. 

The reconstructed zenith-angule and energy distribu- 
tions of data are shown in Fig.l (dots). Eight events were 
reconstructed as upward going cascades (zenith angle 
9 > 90° in the left panel and the distribution in the 
dashed box in right panel). Also the MC-generated (his- 
tograms) and reconstructed (boxes) zenith-angule and 
energy distributions from simulated atmospheric muons 
are shown in Fig.l; 12 upward reconstructed cascade- 
like events are expected. As seen from Fig.l, within sys- 
tematic and statistical uncertainties there is no significant 
excess above the background from atmospheric muons. 
We introduce the following final neutrino signal cuts on 
the cascade energy: E s h > 130 TeV and E s h > 10 TeV for 
downward (40° < 8 <90°) and upward going cascades, 
respectively (see for details lfl6l ). Furthermore, events 
which fulfil selection requirements used in our previous 
analysis 01 are also considered as neutrino candidates. 
With zero observed events and 2.3±1.2 expected back- 
ground events, a 90% confidence level upper limit on 
the number of signal events of n 90 % =2.4 is obtained. 

A model of astrophysical neutrino sources, for which 
the total number of expected events, N mo deh is larger 
than n 90 %, is ruled out at 90% CL. Table I represents 
event rates, detection energy range which contains the 
central 90% of the expected signal, median energy 
of the expected signal, and model rejection factors 
(MRF) n 90 %/N mo dei for models of astrophysical neu- 
trino sources obtained from our search, as well as model 
rejection factors obtained recently by the AMANDA 
collaboration ||2), 0. The model by Stecker [17] labeled 
"S05", represents models for neutrino production in the 
central region of Active Galactic Nuclei. Further shown 
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V :V :V =1:1:1 
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Fig. 2. All-flavor neutrino flux limits and theoretical bounds (see 
text). 



in Table I are models for neutrino production in AGN 
jets: calculations by Protheroe [18] and by Mannheim 
fl9l . which include neutrino production through pp and 
P7 collisions (models "P ^7" and "M pp + pj", respec- 
tively), as well as an evaluation of the maximum flux due 
to a superposition of possible extragalactic sources by 
Mannheim, Protheroe and Rachen |20| (model "MPR") 
and a prediction for the diffuse flux from blazars by 
Semikoz and Sigl ll2Tl "SeSi". As can be seen from Ta- 
ble I the model "P jry" is ruled out with n 90 %/N mo d e i = 
0.5. 

For an E~ 2 behaviour of the neutrino spectrum and 
a flavor ratio v e : : v T = 1 : 1 : 1, the 90% C.L. 
upper limit on the all-flavor neutrino flux obtained with 
the Baikal neutrino telescope NT200 is: 



E 2 <f> < 2.9 x 10" 7 cm-V 



_1 GeV, 



(5) 



for 20 TeV < E„ < 20 PeV. Fig. 2 shows our upper limit 
on the all-flavor E~ 2 diffuse flux, which is a significant 
improvement of the earlier obtained limit Also 
shown are the limits obtained by AMANDA 0, Q 
and Pierre Auger Observatory 11221 . theoretical bounds 
obtained by Berezinsky f23l , by Waxman and Bahcall 
Il24ll . by Mannheim et al.(MPR) J20), as well as the 
atmospheric conventional neutrino fluxes ll25ll . 

IV. Conclusion 

The neutrino telescope NT200 in Lake Baikal is 
taking data since April 1998. Due to high water trans- 
parency and low light scattering, the detection volume 
of NT200 for high energy neutrino events is several 
Megatons and significantly exceeds the instrumented 
volume. This results in a high sensitivity to diffuse 
neutrino fluxes from extraterrestrial sources - more than 
an order of magnitude better than that of underground 



searches and similar to the published limits of the 
AMANDA neutrino telescope. The upper limit obtained 
for a diffuse all-flavor neutrino flux with E~ 2 shape is 
E 2 <S> = 2.9 x 10- 7 cm- 2 s- 1 sr- 1 GeV. 
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